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^ ; Abstract 

O ' 
O . 

■ We derive the two-particle fluctuation correlator in a thermal gas of vr-mesons 

. to the lowest order in an interaction due to a resonance exchange. A dia- 

\ grammatic technique is used. We discuss how this result can be applied to 

ly-^ \ event-by-event fluctuations in heavy-ion collisions, in particular, to search for 
the critical point of QCD. As a practical example, we determine the shape of 

'kj^ ■ the rapidity correlator. 
I> ■ 

O ; I. INTRODUCTION 



o 



The goal of heavy-ion collision experiments is to shed light on thermodynamic properties 
of strongly interacting matter and to determine the phase diagram of Quantum Chromody- 
namics (QCD). Recently, experimental |T|-^ and theoretical |^^T| study of event-by-event 
fluctuations has attracted interest, in particular, because these fluctuations carry informa- 



g_) \ tion about QCD thermodynamics. Recent experimental data from RHIC and CERN SPS 
' show non-trivial pattern of fluctuations which awaits theoretical understanding. 

Thermodynamic description of the final stage of the heavy-ion collision is motivated, to 
^ . a large extent, by a remarkable phenomeno logical success of statistical thermal model in 
^ ! describing the observed particle spectra p2|. As far as the fluctuations are concerned, the 



observed fluctuations are Gaussian to a very good approximation, and are consistent with 
being largely of thermodynamic origin (see for more quantitative analysis). In this paper 
we shall not discuss the relationship between the thermodynamic and observed fluctuations 
in heavy-ion collisions in detail (in particular, we shall not discuss possible importance of 
non-equilibrium effects — see, e.g., [^). At temperatures of order 120 MeV characteristic 



of the flnal stage (freezeout) of a high-energy heavy-ion collision, hot QCD is, to a good ap- 
proximation, a gas of predominantly pions interacting through exchange of resonances. This 
motivates the main goal of this paper: to calculate equilibrium thermodynamic fluctuations 
given the interaction of pions. 

The most interesting application of our results is the search for the critical point of 
QCD. As discussed in PJ^, such a point is characterized by a divergent correlation length, 
i.e., a massless excitation, in the channel with the quantum numbers of the sigma meson 
(scalar isoscalar). The strategy of the search is to scan the QCD phase diagram by varying 
experimental parameters, such as the energy of the collision, looking for the signatures of 
the critical point 0,0] . The fluctuations of pions are sensitive to the presence of light sigma 
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excitation because of the direct coupling of pions and sigma: aim. Thus, measuring the 
pion fluctuations one could determine if the freezeout is occuring near the critical point [0 . 
Another important application of our results is the study of charge fluctuations p!^-p!7 . 



The purpose of the paper is more generic. We wish to derive a formula for the pion 
fluctuations, specifically, for the two-pion correlator, given the effective interactions of pions. 
We shall discuss in detail the effect of the interaction induced by the sigma exchange, both 
because of its application to the search for the critical point and because this interaction 
is the simplest. It is relatively straightforward to generalize the analysis and the results to 
interactions induced by exchanges of other resonances, such as p-meson. 

The goal of this paper is to motivate and derive formulas such as (|16]). As a simple 
example of an application of these results we determine the rapidity correlator (|20|) . Other 
applications and a detailed comparison of these results to experiment are deferred to further 
work. 



II. TWO-PARTICLE CORRELATOR AND FLUCTUATIONS 

Most of the commonly used fluctuation measures, such as Gaussian widths of event- 
by-event distributions of mean transverse momenta, ((Apj-)^), or correlation coefficients 
such as, e.g., (ApxAN) {AN is the multiplicity fluctuation), can be related to the two- 
particle correlator, or the two-particle momentum space density, as emphasized in [0,^ . For 
completeness, and to underline the importance of the two-particle correlator, we shall repeat 
some typical relationships here. An expert may wish to skip this section. 

Let us denote by (...) an average over events. Let us assume that we binned the phase 
space into bins labeled by the values of momenta p (or k). We denote by Up the occupation 
number of a given bin in a given event. Then Aup = Up— (up) is the fluctuation, in a given 
event, of this occupation number around the all-event mean. The two-point fluctuation 
correlator which will be discussed and calculated in this paper is then given by 

(AnpAuk). (1) 

For example, the event-by-event fluctuation of the mean transverse momentum is related 
to the correlator ([^): 

{{Apt)') = -^Y.Y.{^^P^^'^)iPT - W){kT - W), (2) 

' P k 

where pr is the inclusive mean transverse momentum. Similarly, the correlator of pt and 
multiplicity fluctuations (also known as linear correlation coefficient), can be expressed as 

(Ap^AA^) = J^^^(AnpAnO(pT -|¥). (3) 

/ P k 

The multiplicity fluctuation is related to (P by an obvious formula: 

{{ANr) = j:j:{^npAn,). (4) 
p fc 
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It is straightforward to generalize these relationships to correlations of fluctuations of 
momenta or multiplicities of particles of different species, e.g., vr"*" and vr". The most inter- 
esting example is the charge fluctuation, AQ, which can be expressed via a correlator such 
as®: 

((Agf) = EEQV(An^Anf), (5) 

ap i3k 

where is the charge of the particle labeled a. 

The correlator (|I]) also equals {npUk) — {np){np,). Since {upTip,) — {np)6pk = p2{p,k) is 
the two-particle momentum space density, one can see that the absence of correlations, i.e., 
P2{p, k) = p{p)p{k) translates into 

{AnpAnk) = {np)5pk. (6) 

In the literature such trivial fluctuations are termed "statistical". For example, with (^: 
((AA^)^) = (A^), (AptAN) = 0, etc. Corrections to (|^) are commonly termed "dynamical" 
fluctuations. In this paper we shall calculate corrections to (P) due to interactions. 

The relationships of the type (|])-(|) and (|^) demonstrate that the correlator (|l|) is a 
very universal quantity. 



III. CALCULATION OF THE CORRELATOR 

A. Correlator in an ideal gas 

Our goal is to calculate the universal fluctuation correlator which we define as 

i^K^ni), (7) 

where (. . .) denotes the average over a thermodynamic ensemble. An" = rip — (rip) and rip 
is the occupation number for a particle of type a in the momentum mode p. In view of 
experimental applications we shall focus on charged pions, i.e, a = +,— for vr"^ and 7r~. 



In the ideal Bose gas this correlator is given by a well-known formula p4 

{An^p An^)free = 5"^V/p(1 + /p)- (8) 

where 

U = K) = (9) 

is equilibrium Bose-Einstein distribution function. Eq. (§) simply means that occupation 
numbers in two different modes are uncorrelated and in each mode the mean square variance 
is as in the Poisson distribution, (up), times the Bose enhancement factor, (l + rip). Without 
the Bose enhancement factor the correlator is just the trivial correlator (^). The corrections 
due to Bose enhancement, significant in the region of p and k very close to each other, are 
the subject of Hanbury-Brown-Twiss (HBT) interferometry studies and will not be discussed 
here (see, e.g., p5[ for review) 



It is useful to note that /p(l + fp) in the r.h.s. of (|) is the first derivative of f{[3ujp): 

fp^nf3oop) = fp{i+fp). (10) 
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B. Effect of interaction 



We wish to find the effect of the interaction on the correlator (^. Once the interaction is 
turned on we must carefully rethink the meaning of Up. The spectrum is no longer a direct 
superposition of the single-particle spectra, i.e., the energy of a level is not just J^p^p^p- Iii 
other words, the scattering now makes occupation numbers of individual momentum modes 
bad quantum numbers — they are not conserved. However, if the changes in the spectrum 
are controUably small, we may be able to trace the shifted levels to their original position 
in the free gas and assign their quantum numbers correspondingly. We shall see below how 
this qualitative picture works in the lowest non-trivial order in the interaction. 

Such an "adiabatic" definition of the occupation numbers would correspond to a situation 
in which the interaction (i.e., scattering) in a thermalized system is turned off at a point 
in time, the momenta are thus frozen, and the particles are allowed to stream freely into a 
detector surrounding the system. Such a situation can be considered as a model of freezeout 
in heavy-ion collisions. 

Such a picture relies on the weakness of the interaction and, presumably, it cannot 
be pushed beyond the leading order calculation systematically. However, we would like 
to apply it to the gas of rather strongly interacting particles (such as pions). We shall 
argue that some of the contributions of interactions can be absorbed into self-energy, or 
vertex type renormalizations. We then assume, unfortunately without a proof, that the 
resulting effective interactions, to lowest order, would be sufficient to determine the correlator 
to a satisfactory (from experimental point of view) precision, in the same way as a tree 
level effective Lagrangian describes strong interaction amplitudes. Although our motivation 
requires us to make such an assumption, in the remainder of this section we shall simply 
concentrate on performing a formal lowest order calculation, assuming that coupling constant 
can be made controllably small. We shall return to the discussion of applicability of our 
results again. 

C. Interaction and free energy 

In summary, this is the strategy. The energy of a level, descendent from a free level with 
occupation numbers given by a set {up} is i?[{np}] = J^p^pUp + SE[{np}], where 6E is a 
correction which depends, in principle, on all Up. Once we know the spectrum £'[{np}] we 
can determine the joint probability distribution exp{—PE[{np}]) for the sets {up} and thus 
calculate correlators such as (|^). [] 

Rather than pushing forward along these lines, we shall, instead, begin by a more familiar 
calculation of the free energy of the interacting gas using perturbative expansion. We shall 
see that this calculation can be mapped onto the scheme outlined above. 



One should note here that such a calculation must break down if there is mixing between levels 
which are degenerate without interaction (scattering or decay of an unstable resonance is related 
to such a mixing). In this case, 6E becomes singular in the leading order in perturbation theory. 
This fact we shall observe in our calculation. 
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We consider a theory with two charged pions vr"*" and vr of mass m and a scalar field a 
of mass m, with the interaction controlled by a coupling G and given by the Lagrangian 

C^^^ = 2Gan+n-. (11) 

We shall only concentrate on the effect of this simple interaction. It will then be rather 
straightforward to include other fields and their interactions. 





(a) (b) 
FIG. 1. Diagrams for the free energy at order G^. 

To order G^, the correction to the free energy is given by the two diagrams in Fig. |I[ We 
write the corresponding expressions using the mixed time-momentum representation (see, 
e.g., p6|-|2^) for the propagators: 

A(r, p) = (e-'^^" + e-"''(^-")) . (12) 

In this representation the integral over the time r separation of the vertices in these diagrams 
is trivial. The result can be written as: 

F^'^^ = -2G^ / / 7^7^— + 2/p)(l + 2/fc) , (13a) 



1 1 1 



E 



p Jk 2uJp 2u!k 2uJq |^__|_]^| J2i ^i^i 



2 



l + 3i 



(13b) 



where Jp = J (i^p/(27r)^, q = p + k. The first expression, F^'^°'\ is straightforward, while 
for F^'^'^^ we introduced shorthand notations: i = 1,2,3, cui = cUp, uj2 = uj^., tu-^ = tDg, and 
analogous notations for /j's. To elucidate the expression ( |13b| ), we write one term from the 
8 terms in the sum over sets {sj} explicitly. The term 



-4G' 



p Jk 2ujp 2uik. 2l-i^q l^p + CUfc + iiJq 



;i + /p)(i + /fc)(i + /,) 



(14) 



corresponds to {sj} = {1,1,1}. The meaning of this term is simple. It is the ensemble 
average of the second order perturbative correction to the energy, since to leading order 
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5F = {5E) =Y,5EAe-^''^/Z. 

A 



(15) 



The sum cUp + cu^ + in (0) is the energy denominator. The term exhibited in (|T4D comes 
from the virtual intermediate state A' = {n^ + + 1,'n-q + 1}, which has three extra 
particles 7r+, tt" and a, relative to a given state A = {n^,n^,nq}. All terms in ( p.3b|) 
correspond to virtual transitions generated by 2G'cr7r+7r^. There are 8 variants: a state with 
either one more or one fewer a (53 = +1 or —1), and either one more 7r+ or one fewer 7i~ 
[si = ±1), and either one more 7r~ or one fewer 7r+ (s2 = ±1)- The first, simpler, diagram, 
pC^a) ^ describes the contribution of the virtual states with one more/fewer a, and unchanged 
numbers of 7r+ and 7i~ . 

We see that the energy of a given state characterized by a set of occupation numbers 
A = {np,n'^,hq}, receives a correction 6E[{np,n'f^,hq}], which is a function(al) of the 
occupation numbers A in the unperturbed state. Thus, at least up to order G^, we can view 
the energy Ea = Ef + 5Ea of a given state as a function(al) of the occupation numbers in 
this state. This allows us to evaluate thermodynamic averages such as (^. 

D. The correlator 

At zeroth order (i.e., in the free theory), since E^a is linear in np, exp {-l3Ef) factor- 
izes into Gibbs distributions for individual momentum modes, and the correlator is given 
by (H). However, the correction exp(— n^, nq}]) does not factorize, introducing 
correlations between occupation numbers of different modes. 

One can derive simple rules which allow to read off the correlator directly from the ex- 
pression for the free energy. To facilitate this, one can introduce "chemical potentials" for 
relevant occupation numbers, shifting the energy E[{np,n^,hq}] by — Z^a Z]p A'-p'^p- Differ- 
entiating resulting /i-dependent free energy with respect to /i" and /if we obtain correlators 
such as (|^. In the zeroth order this, of course, reproduces (H). In the order G^ the result can 
be represented by diagrams in Fig. These diagrams are obtained from those in Fig. |l|. A 
cross represents differentiation with respect to /i" and amounts to replacing corresponding 
/p or (1 + /p) factors with /^(l + /"), which is the derivative of with respect to [3ujp 
(see (0)). g 



^Although, in equilibrium, and /p are the same, we are using different notations for them 
since they differ when /ip and ^~ are introduced. 
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(a) 



(b) 





(c) (d) 
FIG. 2. Diagrams for the correlator 

There is a class of diagrams in Fig. |^ which affect the correlator (0) in a rather trivial 
way. These are the diagrams (c) and (d), which have both crosses on the same propagator 
line. Such diagrams give contribution proportional to S°'^6pk. These diagrams represent 
O(G^) self-energy type corrections 6ujp to Up multiplied by a factor fp". Since at zeroth 
order the correlator is 5"^5pkf'{(3ujp) (|^), one can see that such diagrams are accounted for 
when we replace Up by Up + 5ujp in f\[3ujp). 

We shall not study the effects of such self-energy contributions on the correlator or on 
the single-particle distributions (np).| Instead, we assume that these corrections have been 
already accounted for and we are using corrected Up. We can view this procedure in the sense 



of the Wilsonian effective theory approach, and consider the Lagrangian with interaction (11 



as the effective Lagrangian. That means no loops (such as self-energy) should be computed 
with it. From a practical point of view, one can imagine that (rip) is measured rather than 
calculated, and the results are expressed in terms of such measured distribution, instead 
of fp- 



^It is useful to note that 0{G'^) corrections to (rip) can be also calculated in the same way. These 
are given by the same diagrams but with one cross only, i.e., proportional to fp. One sees that 
these corrections are accounted for by replacing cOp with its 0(G^)-corrected value in the function 

fiPiOp). 
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With this understanding we can now focus on the most interesting contributions to the 
correlator, given by the diagrams in Fig. |(a,b). Differentiating (|l^) with respect to /ip and 
/if, and discarding self-energy terms proportional to 8pk, we find the order correction to 
the correlator: 



2 fpfk 



+ - 



7' 



al3 



[UJp + Uk) 



+ — 



(16) 



where f'p = fp{l + fp), 7+^ = 7^+ = 1, 7++ = 7 = 0. The first term in the square 
brackets is from the diagram in Fig. 0(a), while the second and the third are from Fig. |^(b). 
A simple way to visualize the a/3 and pk dependence in these three terms is to consider a 
process of forward scattering TTpTrf — >■ VTpTr^. The three terms arise from contributions with 
the sigma propagator in t, s and u channels respectively, as illustrated in Fig. 0. 






t 

FIG. 3. A representation of 3 terms in the correlator (|l6|). The resonance is exchanged in t, s 
or u channel. 



IV. SUMMARY AND DISCUSSION 

In this paper we considered pion gas in thermal equilibrium as an idealization of the 
final (freezeout) stage of a heavy-ion collision. We considered the effect of pion interaction 
on the fluctuation correlator (|I|) or (|^ — an important universal quantity characterizing 
pion fluctuations. Eq. (|TBp gives the lowest order correction to the correlator (|^) due to the 
interaction mediated by the exchange of cr-resonance and given by the Lagrangian (0). Our 
main motivation to consider such an interaction is the search for signatures of the critical 
point on the QCD phase diagram, the point where the sigma mass vanishes. One can see 
that the first term in ([T6|) , the one due to the sigma exchange in t-channel, dominates in the 
limit m — s> 0. This term has been already calculated in [0] using a different method. The 
advantage of the diagrammatic approach of this paper is that it allows one to obtain other 
terms (due to s- and w-channel exchange), subleading in the regime m ^ 0. 

Strictly speaking, we have calculated the correction to the correlator only in the leading 
order of perturbative expansion. Thus, our results are only rigorous in the limit of small 
coupling constant. Our motivation, however, has been the calculation of a two-particle cor- 
relator in hot QCD, which, at temperatures of around 120 MeV, characterizing the freezeout 
stage of the heavy-ion collisions, is primarily a gas of pions, whose interactions are not con- 
troUably small. We observed that, at the leading order, some of the contributions could be 
naturally absorbed into the self-energy renormalizations. We thus conjectured that if one 
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uses an effective (rather than bare) interaction in (0) and effective self-energies, our leading 
order calculation could be sufficient in practice. 

Below we make several comments and outline some questions and problems, related 
to the result (|16D and its comparison to experimental data. More detailed study of these 
questions is left to further work. 



A. Other resonances 



The interaction due to the a meson exchange is simpler than interaction due to other 
resonances, such as p, but generalization is not difficult. For example, given the effective 
p^TT^TT^ interaction in the form |^ 

-CpTTTT = igpTvAp'^T^^d^Tc' + h.c, (17) 

one obtains the contribution of the p-exchange in the t-channel in the form 

(An^ An^),.,_„ = Pg^^ip ■ fc)^^^^^^. (18) 

uJpUJk rrip 

Two notable differences from the similar t-exchange contribution of a in ([T6| ) are worth a 
notice. Both are due to p being a vector particle. First, there is an angular correlation 
between the pions due to the [p ■ k) factor. Second, the sign of the correlation now depends 
on the relative sign of the two pions, due to a minus in 6""^ — 7"^. 

Similar formulas can be derived for other resonance-mediated interactions of pions. It 
should be clear that the contribution of the t-exchange of cr-meson is dominant when m„ 0, 
i.e., near the critical point. It can be further distinguished from other resonance contribu- 
tions, such as p, by its charge-independent character (5"^ -|- 7"^ = 1 for all a/3). Most 
importantly, this contribution will show a non-monotonic behavior if the critical point is 
approached and then passed during the scan of the QCD phase diagram. 



B. Resonance exchange vs decay 

It is important to realize the difference between the result (|TBp and a more common 
approach to the correlations due to resonance decays after freezeout |]7|JT^ . It is easiest to see 
this difference in the charge dependence of the correlator ([T6|) . The term which has the same 
charge dependence as the contribution of the resonance decays is the one due to s-channel 
exchange: a , or p", decays into tt+tt^, therefore the correlator must be proportional to 7"^. 
The contribution of t and u channel exchanges is entirely different. In order to compare to 



^We should remember that the tree-level interaction which determines the pion correlator is the 
effective interaction, which should includes all thermal self-energy and vertex corrections. There- 
fore, the magnitude of the coupling is not necessarily the same as in vacuum (although we do not 
expect it to be significantly different in the case of p-meson), and the Lorentz invariance, which we 
assume for simplicity, is not exact. 
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experiment, however, one has to take into account the fact that the detected pions consist to 
a large extent (more than a half) of the products of resonance decays. The contribution of 
such resonance pions to the correlator should be calculated separately, using the kinematics 
of the resonance decays. Such a calculation has been done numerically, e.g., in for 
fluctuations, where the effect turned out to be small, except for the contribution of the decay 
of light a meson near the critical point, and in |jl2| for charge fluctuations, where the effect 
is of order 20-30%. The emphasis of this paper is on the correlation induced by virtual 
resonance exchange among the direct pions. 



C. s-channel pole and resonance width 



Generalization to heavier resonances, such as p, presents another problem. If we consider 
p and k such that their invariant mass is equal to m, the second term in (|16D becomes 
infinite (the s-channel resonance is on mass shell). The perturbative expansion breaks down.|^ 
However, we know that there is no true pole singularity in this case, since the resonance 
has a finite lifetime. Thus one would expect that careful treatment of higher order terms 
should remove this pole divergence and replace it with a smooth Breit-Wigner-type curve. 
Perhaps, it will also allow to treat resonance exchange and decay within the same formalism 
(cf. previous subsection). 

The s-channel divergence does not arise in the application to light a particle near the 
critical point, since in this case < 2m^. 



D. Longitudinal expansion and rapidity correlator 



In this paper we considered an isotropic pion gas (a fireball) at rest. Before one could 
compare our results to fluctuations observed in heavy-ion collisions it is necessary to take 
into account the effect of expansion. In this paper we shall only consider longitudinal 
expansion. It can be taken into account using Bjorken's boost-invariant hydrodynamic 
model. In this approximation, one considers a superposition of fireballs placed continuously 
over large (ideally infinite) rapidity interval. In order to obtain the correlator in this case 
one needs to take the correlator for a gas at rest (Anp Anf )i.cst7 boost momenta p and k 
with rapidity y: {An^^y Anl^y)^^^^ and integrate over y: 

{Anl An^)iong.exp. = / dy Up^y Ukg^y {An'^y An^^y)^^^^. (19) 

J — CXD 

The factors Up, Uk appear because of the non-uniform density of modes in rapidity space, 
i.e., d^p = ujpd^pj,dyp {yp = aXcinh.{pJujp), pj, = {px,Py))- 



^This is due to the degeneracy or level-crossing as mentioned in an earlier footnote. Unperturbed 
energy levels, different by replacing p with a vr^vr^ pair, are degenerate in this case. 

^However, higher order thermal corrections still induce imaginary part of the a propagator. 
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As a practical example, one can use formula (|T9D to determine the form of rapidity 
correlations. For this purpose one needs to fix rapidities Up and y^, and integrate over 
transverse components and kx'- 

^^"^(2/p,l/fc) = ^^ / / (An^An^),o„g.e.p., (20) 

where /p^ = J d'^prp/{2TT)^. The normalization factor, {dN / dy) / {Jp fp) , is such that 
C{yp,yk) = P2{yp,yk) - p{yp)p{yk), where p = dN/dy and p2 is the density of pairs per 
unit rapidity squared. 

Using expression (|T9D in (|20|) one finds, of course, that correlator C{yp,yk) is a function 
of yp — yk only. The dependence on yp — yk can be found substituting (p!6|) for (An^ A?T,^)rest 



in ([T9|) . The detailed form of the correlator depends on the parameters G, m, /? = 1/T. A 
rather robust feature of C{yp — yk), however, is that it vanishes exponentially with increasing 
IVp ~ Ukl on the scale of 1-2 units. This is because of the fast exponential falloff of the 
occupation number factors fpfk in QTBI ) when either p or k becomes large compared to T. 
As an example, the contribution to the rapidity correlator from the t-channel exchange term 



only, i.e., first term in ([Tq) , is shown in Fig. 



correlator 




1 2 3 4 

FIG. 4. A typical contribution to rapidity correlator, y = \yp — y^l- In this example, T = 120 
MeV, G/rh = 4.5 (see for motivation). Plotted is the correlator C in (^) divided by dN/dy. 

Another interesting consequence of ([l6|) is a specific dependence on multiplicity dN/dy. 
For example, the centrality dependence of the freezeout parameters is known to be rather 
weak. Neglecting it, we could predict, using (|T6|) , that the centrality dependence comes only 
from the factor dN/dy, and is, therefore, linear. 
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